Properties of superconductor - Luttinger liquid hybrid systems 
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In this paper we review some recent results concerning the physics of superconductor - Luttinger 
liquid proximity systems. We discuss both equilibrium (the pair amplitude, Josephson current, and 
the local density of states) and nonequilibrium (the subgap current) properties. 
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I. INTRODUCTION 

The properties of a normal metal (N) which is in good 
electric contact with a superconductor (S) are strongly 
modified, a phenomenon known as the proximity ef- 
fect Jl) . This effect is due to the presence in N of Cooper 
pairs leaking from S, thereby giving rise to a nonvanishing 
local pair amplitude. The microscopic origin for charge 
transfer across the N-S interface is the phenomenon of 
Andreev reflection . The distance over which the pres- 
ence of the superconductor is felt in N is determined by 
the length £jy over which the two electrons forming the 
pair remain correlated. This length decreases with in- 
creasing temperature T. In particular, for clean normal 
metals, £/v = vf/T, where vp is the Fermi velocity. The 
proximity effect manifests itself in various ways. An ex- 
ample is the Josephson effect, occurring in S-N-S sand- 
wiches as long as the thickness of the N-layer does not 
exceed £jv |§- Another example is the local single parti- 
cle density of states (DOS) of the metal, which acquires 
an energy dependence similar to the well-known BCS- 
DOS of a superconductor up to distances of the order of 
£jv away from the N-S interface ||||. 

Due to the recent development of superconductor- 
semiconductor (S-Sc) integration technology, a revived 
interest arose in the properties of clean proximity sys- 
tems. Present-day high-quality Sc heterostructures com- 
bine a number of attractive low-temperature electronic 
properties. The elastic mean free path can be as long as 
20/im. In typical low-density systems, the Fermi wave- 
length \p is of the order of 50nm. These lengths are 
much larger than the corresponding ones in an ordinary 
metal. In addition, the dephasing length L^, over which 
the phase coherence of a single electron is maintained, 
can easily be of the order of 40/im. Thus, quantum bal- 
listic electron propagation dominates in small-scale Sc 
systems. The interplay between phase-coherent electron 
propagation in Sc and macroscopic phase coherence in S 
gives rise to interesting new physics [^^|. 

During the past years, an increasing interest devel- 



oped in the effects of electron-electron interactions on 
the properties of low-density Sc nanostructures. The key 
point is that the Coulomb energy can become comparable 
to the kinetic energy, i.e., the parameter e 2 /vf ~ ^f/^b, 
where e is the electron charge and as is the Bohr ra- 
dius, is no longer small. Moreover, screening becomes 
less effective as the system dimensionality is decreased. 
As a result a non-perturbative, microscopic treatment 
of interactions is required. For one-dimensional (ID) 
systems this can be done in the framework of the Lut- 
tinger model j[o|. In a ID interacting electron system, 
also referred to as a Luttinger liquid (LL), there are 
no fermionic quasiparticle excitations. Instead, the low 
energy excitations of the system consist of independent 
long-wavelength oscillations of the charge and spin den- 
sity which propagate with different velocities. The prop- 
erties of such a system therefore are strikingly different 
from those of a non-interacting ID system. The Lut- 
tinger model is believed to be relevant for a description 
of transport in a number of physical systems like Sc quan- 
tum wires and edge states in the quantum Hall effect [ pd| . 
We should also mention the recent rapid advances in 
controlled fabrication of single- wall nanotubes . Lut- 
tinger liquid behavior is expected in metallic nanotubes 
with two gapless ID modes of excitations E3]. 

At present, relatively little is known about the influ- 
ence of electron-electron interactions on the proximity 
effect in clean mesoscopic N-S systems. In view of the 
above one may conclude that S-Sc heterostructures are 
good candidates to study such effects. Of particular in- 
terest are ID quantum wires, connected to a supercon- 
ductor (S-LL systems). Various techniques are available 
nowadays to confine electrons in a semiconductor to a 
long and narrow channel. In combination with state-of- 
the-art S-Sc integration technology, this will make sys- 
tematic studies of interaction effects in S-LL systems fea- 
sible in the near future. In this paper we review various 
properties of S-LL proximity systems. The paper is or- 
ganized as follows: In Section p we define the spin-1/2 
Luttinger model. In Section [II, we discuss two possi- 
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ble ways to couple the LL to a superconductor: via a 
highly transmissive and via a tunneling interface. The 
remainder of the review will be divided into two parts. 
In Sections [v], and VI we discuss equilibrium prop- 
erties of S-LL heterostructures: the pair amplitude, the 
Josephson effect, and the local density of states. In Sec- 



tions VII and VIII we will focus on transport phenomena 



and calculate the subgap current through a S-LL inter- 
face. 



II. THE SPIN-1/2 LUTTINGER LIQUID 

In this Section we set the notation which is needed 
in the rest of the paper. For a detailed discussion of 
Luttinger liquids we refer to existing reviews on the 
topic jlfl]. The long- wavelength Hamiltonian of a ID 
interacting electron system of length L can be expressed 
as that of a harmonic fluid for the charge (j — p) and 
spin (j = <y) degrees of freedom 



Ht 



dx 
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The parameters gj are related to the interaction strength 
(gj = 2 for non- interacting electrons), and Vj = 2vp/gj 
are the velocities of spin and charge excitations. The 
commutation relation between the Bose fields 9(x) and 
4>(x) for each spin sector (spin up s = +1 and down 
s = -1) is [<p s (x),6 s t(x')] = (in/2) sign(x / - x)5 a a r, 



where 



s0 CT) and 6 a = P + s9„ 



The parameters gj can be determined once one defines 
an appropriate microscopic Hamiltonian. For a Sc quan- 
tum wire with spin-independent interactions, one may 
take g a = 2 and g p = 2/y/l + 2Vq/itvf, where Vq is 
the zero-momentum Fourier component of the interac- 
tion potential |l4|] . Spin and charge excitations propagate 
with different velocities (spin-charge separation) . In gen- 
eral there will be additional nonlinear terms appearing 
in Eq. (Q) due to backscattering or Umklapp processes. 
Usually, in quantum wires away from half filling and in 
the ballistic regime these terms can be ignored; therefore 
we will not discus them in the following. The electron 
field operator ^ is expressed in terms of the boson field 
in terms of the spin and charge degrees of freedom 



*.(a:) 



<5=± 



(2) 



where kp is the Fermi wave vector and pa = kp/n the 
electron density. 



III. SUPERCONDUCTOR-LUTTINGER LIQUID 
INTERFACES 

A superconductor can contact a quantum wire at the 
end (edge contact) or at some internal point or segment 



of the wire (lateral contact). Both edge and lateral con- 
tacting have been already implemented in experiments 
with S-Sc structures The quality of S-LL interfaces 
can be characterized by the transparency of the barrier 
at the interface. We will consider below edge contacts 
with both high and poor transparency of the barriers as 
well as poorly transmitting lateral point contacts. 

Perfectly transmitting interfaces - Maslov et al. (ilf | 
and Takane and Koyama [l^ recently developed a 
bosonization scheme to treat clean S-LL interfaces. Fol- 
lowing Ref. we consider the case of two supercon- 
ductors, kept at a phase difference %, and adiabatically 
connected to a quantum wire of length L (the results 
for a single interface can be obtained by taking the limit 
L — > oo). The bosonization scheme can be carried out 
once the boundary conditions for the right and left mov- 
ing fields at the S-LL interface are determined. 

The mode expansions for the fields 9j and <pj (j — p,a), 
such that the Fermi operators (Q) satisfy the proper 
boundary conditions at the interface, are 




V 9>0 



(3) 
(4) 




<j) a (x) =f a + J— ^2lq{bt,q + K,q ) 
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q>0 



q>0 



where q — irn/L 7+ = (i:/2qL) 1 / 2 cos(qx), and 7 ? = 

i(-K /2qL) 1 / 2 sijx(qx). The system is characterized by the 
two topological numbers J p and M a ; the topological 
numbers N s = (M a ± J p )/2 must be odd integers Jlg[ . 

Substituting the mode expansions @ - (§) into the 
phase Hamiltonian ([!]), we obtain a Hamiltonian of the 
form 



7T 

H =4L 
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j=p,a q>0 
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Takane and Koyama |17|] extended this bosonization 
scheme to include the energy dependence of the phase 
shift related to Andreev reflection, thereby showing that 
it may be important if a sufficiently strong potential bar- 
rier is placed at the S-LL interface. In the following we 
will neglect this energy dependence. 

Poorly transmitting interfaces - If the quantum wire is 
weakly coupled to a superconductor at x = 0, the sys- 
tem can be treated in terms of the tunnel Hamiltonian 
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formalism. The Hamiltonian of the whole system p8| , 
H = Hs + Hl + Ht, contains the BCS-Hamiltonian Hs 
describing the bulk superconductor, the Hamiltonian Hl 
of the quantum wire, and the tunnel Hamiltonian, 



H T = Y, t ^sA x = °)^lAx = 0) + (h.c. 



(8) 



The tunnel matrix elements to can be related to 
the tunnel conductances Gt of the junctions, Gt = 
47re 2 A L (0)A s (0)|t | 2 , where A L (0) = l/irv F , and N s (0) 
is the normal state density of states at the Fermi level of 
the superconductor. 

If two superconductors are attached to the ends of 
a quantum wire one should use the recently developed 
bosonization technique for finite ID systems with open 
boundaries jl9| . The presence of lateral tunnel contacts 
does not impose any additional boundary conditions on 
the Fermi fields. We will neglect possible inhomogencities 
of the quantum wire near the contacts, which might be 
the source of electron backscattering. For a discussion 
of possible modifications due to the presence of a barrier 
potential, see Ref. p0[. 



IV. PAIR AMPLITUDE 

In this section we will evaluate the pair amplitude in- 
duced into a Luttingcr liquid which is connected to a su- 
perconductor at x — 0. We calculate the pair amplitude 
at a distance x from the contact. The pair amplitude is 
defined as the anomalous time-ordered expectation value 
H(a;, t) = —(T T ipf(x,T + )il>i(x,T)), where r + tends to r 
from above. Throughout this section, we will be inter- 
ested in the smooth spatial variation of 3 (we ignore 
spatial variations on the scale of Xf)- We will evalu- 
ate 3 both for a highly transmissive and for a tunneling 
interface. 

Perfectly transmitting interfaces - The pair ampli- 
tude for an adiabatic interface between S and LL 
has been obtained by Maslov et al. 3(x) = 

—2p (exp2i(f)p}(cos29 a ). Since Maslov's bosonization 
procedure has been developed for Andreev scattering at 
energies much smaller than the superconducting gap A, 
we shall evaluate 3 at large distances from the inter- 
face, i.e., with a short wavelength cut-off a such that 
i » a ~ e uj?/A. Performing the averages with 
respect to (M), we obtain 



3(x) = -2p Q 



ira/ '/3vf 
sinh 2ttx//3vf 



1/2 r 



ira/j3vp 



sinh 2irx/[3vp 



(9) 



where (3 — l/T. The result consists of a product of a 
spin (first term in brackets) and a charge contribution 
(last term in brackets); only the latter is sensitive to in- 
teractions. 



At zero temperature, the pair amplitude for non- 
interacting electrons decays slowly away from the junc- 
tion for: Et=o(x) ~ 1/x. Nothing prevents the supercon- 
ducting correlations from being present arbitrarily deep 
in the quantum wire. In the presence of repulsive inter- 
actions, the two electrons can be scattered out of their 
time-reversed state and therefore 3 decays faster with 
increasing distance x from the interface, 



a 
2x 



1/2+1/9, 



(10) 



At finite temperatures, the coherence length in the LL 
becomes finite, and, correspondingly, the pair amplitude 
should decay faster. This can be seen by calculating the 
ratio St/3t = o. At low temperatures T <C vp/x, this 
ratio behaves as: 



5 T (x) ^ _ 1 / nxT Y 
S T=0 (x) ~ 6 V v F J 



(11) 



independent of the interaction strength. At higher tem- 
peratures T vf/x, the suppression will be exponential: 



3 T =o(^) 



,(-2ttxT/v f ) 



(12) 



Here, interactions determine the pre-exponential temper- 
ature dependence. 

Poorly transmitting interfaces - In order to evaluate 
the pair amplitude for a tunnel interface between S and 
LL, we use standard perturbation theory in the tunnel 
Hamiltonian Ht, Eq. (||). In the tunneling regime the 
usual bosonization scheme can be used, with short wave- 
length cut-off a ~ Xf ■ The lowest order non- vanishing 
contribution to the pair potential in the LL is given by 
the following expression 



E(x,t) = % J dT 1 dT 2 Tl + ,-,- t+ (x,T + ]x, t;0, n;0, r 2 ) 



x F +! _(0,-r 2 ;0,Ti). 



(13) 



Here F_|_,_(0, r 2 ; 0, t\) is the usual anomalous Green's 
function for the superconductor, and we intro- 
duced the four-point correlator Hi, 2, 3,4(1, 2, 3, 4) = 
(T T ^ s i(l)^) S 2(2)^i] 3 (3)^>J 4 (4)). The relevant process con- 
sists in the tunneling of two electrons from S into the LL. 
This process is of second order in the tunneling. First, an 
electron tunnels into the LL, leaving behind a quasipar- 
ticle excitation in S. Then, the second electron tunnels 
from S into LL, annihilating the quasiparticle in S. Thus, 
two time scales play a role in this tunneling process: the 
lifetime |n — ra| ~ 1/A of the intermediate state with one 
quasiparticle in S, and the time |r — Ti|, |t — r 2 | ~ |x|/i>f 
during which each electron propagates in the LL to the 
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position x. The behavior of the pair amplitude will de- 
pend on the relative magnitude \x\A./vf of these time 
scales for propagation and tunneling. 

Pair amplitude close to the junction - For distances 
close to the junction, £g ^> \x\ ^ a ~ Xp, the time 
needed for the pair to traverse the LL is negligible com- 
pared to the quasiparticle lifetime 1/A. Hence, the 
main contribution to 3 stems from electrons that prop- 
agate fast and independently through the LL. In this 
case, we can approximate the four-point correlator by 
the product of two Green's functions II ~ G--(x,t — 
ti)G + , + (x,t + - t 2 ) with G s>s >(x,t) ~ 5 S:S >G(x)5(t). 
The function G(x) is obtained by integrating the sin- 
gle particle Green's function G gs r(x, r) over imaginary 
time (see Appendix A). Within this approximation we 
find E(x,t) ~ ^G 2 (x)F + _(0,t+;0,t). The anomalous 
Green's function can be written as F + ^(0,t + ;0,t) — 
iVs(0)Asinh -1 (cj£)/A), where the Debye- frequency ujd 
is a high energy cut-off in S. We thus arrive at 



E(x) 



where 





"5p" 


1/Sp+9p/4 






. 2 - 




j3v F _ 



(14) 



3° = 2p 



G T oA 

(4e 2 /7r) vp 



sinh 1 (u> D /A). 



(15) 



Expressions for the function FG(x,g p ) are given in Ap- 
pendix A. 

In the noninteracting case g p = 2 the pair ampli- 
tude becomes space-independent at zero temperature, 
Eg p= 2(x) — ► 3°. This reflects the fact the time- reversed 
electrons from S reach points in the LL within a distance 
~ £s from the junction instantaneously. In the presence 
of repulsive interactions, propagation through the LL be- 
comes faster, v p > Vp, and one naively would expect the 
pair amplitude to remain space independent. However, 
for g p < 2 we find at zero temperature, 



c0 



S T=0 (x) = 3' 



l/9p+ff P /4-l 



Fh(9 P ) 



(16) 



Due to repulsive interactions, the pair amplitude be- 
comes x— dependent: it decays algebraically away from 
the junction. This suppression of the pair amplitude is 
related to the fact that the effective tunneling amplitude 
of Cooper pairs is renormalized at low energies by the 
interactions, similar to the tunneling amplitude of single 
electrons into a LL [jllj . 

At finite temperatures, we therefore expect a competi- 
tion between two effects: (i) enhancement of the tunnel- 
ing due to thermal fluctuations leading to an enhance- 
ment of the pair amplitude as a function of temperature; 
(ii) decreasing coherence length, leading to a suppression 
of the pair amplitude with temperature. As a result of 



the above competition, the pair amplitude shows a max- 
imum as a function of temperature if g p < 2. Analytical 
results can be obtained in the case of weak interactions 
(1 — g p /2)irx/ f3vp <C 1. At low temperatures (3vf ^> x 



^T=0 



l + (2- 5p ) — 



2 r(i/ + l/2) / x 

(3v F 



Pvp 0F r(i/ + i 



(17) 



with v = 1/4<7 P + 5 P /16 + 1/4. Note that the first term, 
describing the thermal enhancement of the pair ampli- 
tude, vanishes in the noninteracting case g p = 2, as it 
should. 

Pair amplitude away from the junction - If we are in- 
terested in the pair amplitude at large distances from the 
interface, \x\ 3> £s, we can neglect the time 1/A spent 
in the virtual state in comparison with the long time 
x/vp needed to traverse the LL. We thus approximate 
the anomalous Green's function in S with the help of a 
5-function in time as 



F +f _(0, n; 0, r 2 ) » irNsiOMn - r 2 ). 



(18) 



This enables us to perform one integration over imagi- 
nary time in Eq. (|13| ) , and obtain 



3(x) 



irN s (0)t 2 J cItUc{x,t). 



(19) 



The function IIc(3;,t) is defined in Appendix B. 

In the noninteracting case g p — 2, the propagator He 
can be integrated analytically. In the zero temperature 
limit [3 — > co, the pair amplitude is given by 



i{x) = 2p 



Gt ot 



(4e 2 /7r) x 



(20) 



We see that 5 decays slowly away from the junction. 
With increasing temperature, the distance vp/T over 
which two electrons maintain their relative phase coher- 
ence decreases, and in the limit (3vp/\x\ -C 1 the pair 
amplitude decays exponentially with x: 



-g„=2 0) - 4 Po 



Gt 2irct 
(4e 2 /7r) (3v F 



exp{-27Tx/j3v F ). (21) 



We now turn to the case of repulsively interacting elec- 
trons, g p < 2. In the zero temperature limit we obtain 



S(x) = 2p c 



Gt 



(4e 2 /7r) 



2/9p 



Fc(g P ), 



(22) 



where Fc{g p ) is given in Appendix B. We find a power 
law decay of 3 with x; compared to the noninteract- 
ing case, Eq. (p0|), the pair amplitude decays faster 
away from the junction. In the high temperature limit 
f3vp/\x\ <C 1, the decay becomes exponential, 
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E(x) ~ 4p 



(4e 2 /Tr) 



(3v p 



2/Sp 



exp(-2nx/Pv F ). (23) 



Note, however, the algebraic temperature dependence of 
S through the pre-exponential factor. 



V. JOSEPHSON EFFECT 

The Josephson effect is an observable consequence of 
the penetration of the pair amplitude into LL |l8|,[l6|,^0| . 
The calculation of the Josephson current through S-LL-S 
system is analogous to the one for the pair amplitude (see 
Section IV). For perfectly transmitting interfaces Jl6| the 
critical current is equal to its value for non-interacting 
electrons independent of the actual interaction in LL. 
This result stems from the fact that features of LL are ob- 
servable in transport experiments only in the presence of 
backscattering, both in normal fll|| and superconducting 
systems. 

The case of poorly transmitting interfaces has been 
treated in Refs. Jl8| within the tunnel Hamiltonian for- 
malism. An infinite quantum wire coupled to two super- 
conductors by tunnel junctions at a distance x was con- 
sidered. In the lowest order in Ht (||), the DC Josephson 
current Ij = I c sin x shows standard dependence on the 
phase difference x between the superconductors. The 
critical current I r . can be estimated as, 



Gt\G 



T2 



(2e 2 /7r) 



(24) 



where the characteristic energy scale 6 and the scaling 
factor Qc characterizing the decay of the Cooper pair den- 
sity in the LL are given by S — A, 3 C = (a/x) 9p / 4+1 / £, ' ,_:L , 
for x -C £g, and 5 = vp/x, 3 C = (a/x) 2 ' 9 <'~ 1 , for x » £s 
(cf. Eqs. (0), (gi)). Similarly to the pair amplitude, the 
critical current shows maximum as a function of temper- 
ature for both x -C £s an d x 3> £<?. At high temperatures 
T vp /x the Josephson current is suppressed exponen- 
tially (cf. Eq. (|2^)). Numerical estimates Jlq] show that 
for typical experimental parameters the critical current 
( p4| ) is in a few nanoamp range. 

If a DC voltage eV <C 2A is applied between the su- 
perconductors, the Josephson phase difference \ becomes 
time-dependent, \ = loj = 2eV , and the Josephson cur- 
rent oscillates with the frequency loj (AC- Josephson ef- 
fect). For non- interacting electrons the Josephson cur- 
rent I j it) oc I ac sin (u)jt — Xo) acquires additional phase 
shift xo = eVx/vF due to the propagation of electrons 
between the contacts. In the interacting case also the 
amplitude I ac of AC current becomes voltage depen- 
dent. Namely, for moderate interaction the amplitude 
Iac(V) shows pronounced oscillations with the period 
eVb = 2tt/[x/v p — xjv a \ corresponding to 27r difference 
between the phases of charge p and spin a excitations 
pSJ. Therefore, the AC Josephson effect can be used as 



a tool for the observation of the spin-charge separation 
in LL. 

A different geometry - finite quantum wire of length 
x S> £s contacting the superconductors at the ends - 
was considered in Refs. The result by Maslov 

et. al. |Q (3 C = (a/x) 2 < 9p ~ 1 ) was revised by Takane 
[p0|| , who found stronger decay of the critical current, 
3 C = (a/x) 2 ^ 2 ^ 9l '^ 1 \ due to a pinning of spin fluctua- 
tions by superconductors. 



VI. DENSITY OF STATES 

The local DOS is defined through the retarded one- 
electron Green's function of the LL Gn(x,x';t) = 
-i{{^(x,t),^(x',0)^})e(t), as 



1 f°° 

N(x,co) = Im dte lult G R {x,x]t) 

n J-oo 



(25) 



In the case of an infinite LL it is straightforward to com- 
pute this quantity and get fo| N(u) - w (9p+ 4 /fp- 4 )/ 8 . 
Contrary to Fermi liquids, whose quasiparticle residue is 
finite, Luttinger liquids have a density of states which 
vanishes at the Fermi energy as a power law, both for 
repulsive (g p < 2) and attractive (g p > 2) interactions. 
In the non-interacting case (g p = 2) the density of states 
is constant, as for a Fermi liquid pp| , pT| . 

What are the modifications of the local DOS due to 
proximity effect? Here we discuss the space and fre- 
quency dependence of the DOS of a LL contacted at 
x = with a superconductor p^j, which corresponds to 
the limit L — > oo in the mode expansion given by Eqs. (j^) 
- (^|) . In this case only the non-zero modes contribute to 
the local DOS. The correlation function (ip^(x, t)tp(x, 0)) 
can be evaluated using the boson representation (the cor- 
relator is not translationally invariant due to the presence 
of the interface). At small energies the DOS behaves as 



(26) 



The exponent of the DOS is negative (g p < 2) and, hence, 
there is a strong enhancement at low energies whereas in 
the absence of S the LL would show a vanishing DOS at 
the Fermi energy. The enhancement of the DOS occurs 
regardless of the distance x from the interface. The scale 
of the enhacement is set by a space dependent high fre- 
quency cutoff lu c ~ Vp/x. The induced pair amplitude 
in the LL, which is characteristic of the presence of the 
superconductor, decays as a power Jlq] of the distance x 
(see Section [Tv] ) . This profound difference in the space 
dependence demonstrates that the DOS provides differ- 
ent information compared to the proximity effect. The 
reason why the DOS does not approach the well-known 
behaviour of an Luttinger liquid far from the supercon- 
ducting contact is in part related to the fact that we are 
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considering a clean wire. In this case the states in the LL 
are extended and the DOS enhancement does not depend 
on x. 

So far we discussed only the case in which the inter- 
face between the superconductor and the Luttinger liq- 
uid has a high transparency. Let us shortly comment on 
the opposite limit, in which the Luttinger liquid is con- 
nected to the superconductor by a tunnel junction. At 
low energies, we find for the DOS close to the junction 
N S -LL ~ w (9p/2-i)+(i/2 Sp -W8). Although the exponent 
is different from the one appearing in Eq. (|26|), the DOS 
is clearly enhanced. Moreover, also in this case the en- 
hancement is found regardless of the distance from the 
junction. 



VII. TWO-ELECTRON TUNNELING INTO A 
LUTTINGER LIQUID 

In this Section, we will calculate the subgap conduc- 
tance for a LL, connected to S via a lateral tunnel junc- 
tion |2^] . The tunnel current is calculated in the standard 
way as 



J(t) 



i(N L (t)) = -ie([H T (t),N L (t)}). 



(27) 



Here, Nl = J2 S I dxtp\ s (x)ipL,s(x); the time-dependent 
tunnel Hamiltonian is given by 



H T (t) = eieyti o< s (0, t)V>s,.(0, t) + h.c. 



(28) 



where eV is the applied bias voltage between S and LL. 
The time dependent field operators are defined as 

where [ii is the chemical potential. We then perform 
an expansion in the tunneling Hamiltonian (which is 
switched on adiabatically at t = — oo), and obtain 



I(t) = -ie / dtxdt 2 dHG R {t,ti,t 2 M), (29) 



where we introduced the retarded Green's function 

G R {t,t 1 M,h)=iO{t-t 1 )6{t 1 -h)6{t 2 -h) 

x {[[[\H T (t),N L (t)},H T {t{)}, H T (t 2 )}, H T (t 3 )}). (30) 

Using imaginary time techniques, Gr can be expressed 
as a product of time-ordered correlation functions for 
LL and S. As long as the relevant energies (eV, T) are 
small compared to the gap of the superconductor, the 
time-dependence of the anomalous correlations in the su- 
perconductor can be approximated with the help of S- 
functions, see Eq. (|l8|). In the zero temperature limit we 
get 



/ = 24e 2 (p a) 2 V 



(g p /2)^(aeV/v F ) a ^ 
T( ai ) 

{g p /2) a HaeV/v F )^- 1 



T(a 2 + 1) 



(31) 



where a\ = 2/g p +g p /2 and a 2 — 2/g p Note that the cur- 
rent depends on the applied bias in a power law fashion, 
which is common for transport through an interacting 
ID system jyj. For noninteracting electrons I ~ V as 
one expects. For repulsive interactions g p < 2, and the 
dominant contribution to current at low voltages reads 

/ rv, V 2 / 9 " . 

At finite temperatures we have 



G — -G t Rk 



ci 



T 
Ef~ 



C2 



T 
Ef~ 



(32) 



Here, E F is the Fermi energy; c\,c 2 are dimensionlcss 
constants, which can be determined numerically (some 
examples are c\ — c 2 — 1.0 for g p — 2.0, c\ = 1.549, 
c 2 = 2.467 for g p = 1.0, and ci = 1.153 and c 2 = 0.765 
for g p = 3.0). 

In the case of a chiral LL the subgap conductance has 
been studied by Fisher 123]. 



VIII. RENORMALIZATION GROUP FOR THE 
SUBGAP CONDUCTANCE 



In this last section we discuss the subgap current by 
means of a poor man's renormalization group ||l7| , p4| , ^5f 
previously developed for the case of an impurity in a 
Luttinger liquid [ji"l|| . To be specific we consider a semi- 
infinite normal metal for x < and a superconductor for 
x > 0. At x — there is an insulating barrier which 
is modeled by a delta-like potential U(x) = UoS(x). In 
the absence of electron-electron interaction in the normal 
metal, the scattering states, due to the presence of the 
interface with the superconductor, are 



l 

\/ 2 TT 



4>k (x) 



ikx 



-ir a e 



ir a e 



roe 

ik* X 



ikx 



' X 



x < 
x < 0. 



(33) 



The states (f>k and (f>k* correspond to incoming particles 
and holes; r and — ir a correspond to the amplitudes for 
normal and Andreev scattering. For the sake of simplic- 
ity, we make the Andreev approximation: in the absence 
of a potential barrier (i.e., Uq = 0) at the interface, one 
has ro = and r a = 1. For our present case of a delta-like 
potential 
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ro = 



2zz* + i(z + z*) 
(1 + iz)(l - iz*) + zz* ' 
1 

(1 + iz)(l — iz*) + zz* ' 



(34) 
(35) 



where z — mUo/k; furthermore k = ^J2m(n + iui n ), and 
uj n = 2-k(ji + 1/2)T. We work with Matsubara frequen- 
cies. To get the proper retarded Green function eventu- 
ally iuj n —> lu + iS as usual. The effect of the interaction 
is analyzed first in the Born approximation. Consider 
for simplicity the Hartree contribution in first order in 
perturbation theory. The scattering states of Eq. j3^ ) 
generate an oscillating electron density, which in turn 
generates a potential Vh (x) . The correction to the wave 
function can be written in the form 



Slpk,a( X ) 



dyG kiaf 3(x,y)V H {y)^k,i3(y), (36) 



where a and (3 take the values 1 and 2 corresponding 
to the particle and hole degree of freedom, respectively. 
The matrix elements of the Green's function take up the 
values 

^ e ik\x-y\ _|_ r ^ e -ik(x+y) 



ir a e 



—i(kx — lz* y) 



Gk,n(x,y) = 
Gk,i2(x,y) = 
G ki21 (x,y) = -f ; 



(37) 



G k ,2a(x>v) = ~T 



^ e -ik*\x-y\ _|_ r * e ik"(x+y)^ 



The correction to the Andreev scattering coefficient is 
obtained by setting a — 2 and 0=1, and reads 



dy [ G k>2 i (x, y) V H (y) <pk,i (y) 

+ GkMx lV )V H {y)cj)ka]- (38) 



By inserting the expression for the Green's function 
Gk,a,p{x,y) and the non-interacting scattering states 



4>k,f}{x) one gets 5r a = 8r^ + 5r^a\ with 

1-0 

(39) 
(40) 



5rW 



ik 



dyV H (y)e 2tk 'y. 



By expanding the Hartree potential in a Fourier series, 
one obtains an expression for the scattering coefficients 



6rP , 
where 



mr a r f°° dq 1 

■— L^^Pm' (41) 

i^)^ (42) 



Sn(q) = / dye- iqv Sn{x). 



(43) 



Notice that the convergence of the integral is automati- 
cally controlled by the imaginary parts of the momenta 
k and k* . The density Sn(x) of the electron is evaluated 
in terms of the non-interacting scattering states, and can 
be written as 



n(x) = no + 8n(x) 



= -l[Tj2 G ^( x ^)+ T J2 G ^ x ' x n ■ ( 44 ) 

For the Fourier transform Sn(q) we then get 



, . ! f ir ( 2k F + q + i5 

\ ztt \ q + 10 

if* , fq — %k F + iS 

- In 1 

2tt 



Hence 



q + i5 



(45) 



, , f2k F + q + i6\ fq~2k F + i5 

x r m I —7 | - r n In 



\ q + i5 



47rfc J_ 0o 2n ' q - 2k 

o - 



q + iS 



and 



^ ( 2) _ imr a r% f°° dq T ,, , 1 



8r?> = 



-V(q)- 



4?rfc 7-oo 27T q + 2k* 
' 1 q + tS 



. . .2k F + q + i5\ . fq-2k F +iS 
x r ta| —7 I - 77, In 



q + id 

Using the Cauchy theorem, we find 



Sr(1) = ™ / / 2k F + 2k + iS 

° 4?rfc V V 2fc + z<5 



-rlV{2k) In 



/2fe - 2k F +id 
V 2fc + i5 



and 



-2fc* + j<5 
-2fc* - 2fc F +«<5 N 



47rfc 

-rSV(-2A;*) In . 
Keeping only the divergent terms, we are finally left with 



Sr a = r -f^V(2k F )ln(^ 
2irv \ e 



(46) 



As it is usual in ID systems, perturbative corrections 
are logarithmically divergent, signalling an instability of 
the Fermi liquid ground state and the emergence of the 
Luttinger liquid behavior. One can then sum the lead- 
ing logarithmic singularity by deriving and solving the 
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renormalization group equation. Here we derive an equa- 
tion for the Andreev reflection coefficient R a = \r a \ 2 - 
The subgap conductance is obtained via the formula 
G = Aire 2 R a , where the condition R a + R = 1 has been 
used. By defining a = ^j- 2 *^) one obtains the equation 



dR a 



2R a (l-Ra) 



where 77 = InA/e with initial condition 77 = (e 
The solution of Eq. (47) gives 



Ra(e) 



i?a(A/e) 5 



R + R a (A/e) 



(47) 
= A). 

(48) 



Here R a and i?o are the initial values of the scattering 
coefficients. Our analysis has neglected the cutoff de- 
pendence of the interaction coupling, which a more com- 
plete calculation should take into account, as well as the 
Fock terms. In Ref. p0| t his calculation has been car- 
ried out and the result (|48j) remains unchanged provided 



APPENDIX B: TWO- PARTICLE CORRELATION 
FUNCTION 



Next we consider the two-particle correlation function 
Hc{x, r) = n+._._.+ (x, t; x, r; 0, 0; 0, 0), which describes 
particle-particle propagation of a spin singlet in a LL over 
a distance x during a time r. Integrated over imaginary 
time, this quantitiy can be written as 

p/2 

/2 {Ot\ 2 /9p 
dTll c (x,T) = p — (- F c (x,g p ). 
VF 

-13/2 



At zero temperature we find Fc(x,g p ) — + Fc(g p ), with 



(Bl) 



28 with 8 = ^ 

r 1 27TV 



OO 

f dz 


I 


1/2 


I 


l/9p 


I 7T 


l + z 2 _ 




[l + (2z/g p )\ 





APPENDIX A: ONE-PARTICLE CORRELATION 
FUNCTION 

In this Appendix, we present some explicit expressions 
for the integrated single particle Green's function of an 

0/2 

infinitely long LL, defined as G(x) — J dTG s . s (x,T), 

-13/2 

where G s ^(x,t) = -{T r ip s {x, r)^,(0, 0)). The latter 
quantitity is readily evaluated pXj ], using the Bose rep- 
resentation (^|). Upon integration over r one obtains in 
the long wavelength limit 

2 , a \ l/2g p +s p /8-l/2 

G(x) = — sm(k F x) - Fg{x,9 p ). (Al) 

vp \ x / 

In the zero temperature limit, Fq is ^-independent, 



F a{x,g P )^ / —sin 

7T 



arctan(l/z) + arctan(g p /2z) 



1 



1 + z" 



— OO 

1/4 



f 



l/4ff P +s P /16 



F G (g P ). (A2) 



l + (2z/g p )\ 

Note that Fq(2) = 1. For noninteracting electrons, 
FG(x,g p ) can be found explicitly, 

2 



Fg(x, 2) = — arctan[l/ svahfax/ ^Vf)]- 



(A3) 



In the case of weak repulsive interactions, (2 — 
g p )Trx/28vF <C 1 and low temperatures itx/2[3vf <C 1, 
the function Fa(x 1 g p ) can be approximated as 



Fg{x,9 p ) ~ F G (g P ) 



2.x: 



9 P vx 



(A4) 



Bvf 2 Bvf 
Here we dropped terms 0[(2 — g p )(nx / Bvf) 2 + (2 — g p ) 2 



Note that Fc(2) = 1. In the noninteracting case 
Fc(x, g p ) can be calculated analytically also at finite tem- 
peratures, 



F C (x,2) 



(2irx//3v F ) 



^cosh 2 {2irx/8v F ) - I 



(B2) 



In the limit of high temperatures, 2-kx/{8vf) S> 1, Fc 
can be found for arbitrary interaction strength, 



F c (x,g p) )n2 



/2ttx\ 2/9p 

f^J exp(-2*x/pv F ). (B3) 



In the case of weak repulsive interactions and low tem- 



peratures an expansion similar to the one leading to ( A4) 
enables us to approximate 



Fc(x,g P ) i2/^\ 2 / ttx 

Fc(9p) ~ 3\0vf) +( 9p) \Pvf 



2/g P 



(B4) 
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